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^ , ABSTRACT 

^vq ■ We present initial results of a new, systematic search for massive star-formation in 

the host galaxies of the most luminous and probably most massive zj>4 radio-quiet 
quasars (M B < -27.5; jaL„(1450A) > 1O 13 L ). The survey, undertaken at 850/im 
J> ' using SCUBA at the James Clerk Maxwell Telescope (JCMT), has a target sensitivity 

limit of 3(7 ^lOmJy, set to identify sources suitable for detailed follow-up e.g. contin- 
uum mapping and molecular line diagnostics. A total of 38 z<;4 radio-quiet quasars 
have been observed at 850/xm, of which 8 were detected (>3ct) with Ss50p m iS lOmJy 
(submillimetre-loud) . The new detections almost triple the number of optically se- 
lected, submillimetre-loud z^4 radio-quiet quasars known to date. We include a de- 
tailed description of how our quasar sample is defined in terms of radio and optical 
properties. As a by-product of our selection procedure we have identified 17 radio-loud 
' quasars with zJb4. 

There is no strong evidence for trends in either detectability or 850/im flux with 
Q ■ absolute magnitude, Mb- We find that the weighted mean flux of the undetected 

| sources is 2.0 ± 0.6mJy, consistent with an earlier estimate of ~ 3mJy based on more 

c/3 , sensitive observations of a sample z^4 radio-quiet quasars (McMahon et ai, 1999). This 

corresponds to an inferred starformation rate of ~1000M Q yr _1 , similar to Arp220. 
The typical starformation timcscalc for the submillimetre-bright sources is ~ 1 Gyr, 
10 times longer than the typical accretion-driven e-folding timescale of ~ 5 x 10 7 
' years. Our 850/zm detection of the z = 4.4 quasar PSS J1048+4407 when analysed 

\ in conjunction with 1.2mm single-dish and interferometric observations suggests that 

this source is resolved on angular scales of l-2"(6-12 kpc). In addition, we present a 
new optical spectrum of this source, identifying it as a broad absorption line (BAL) 
quasar. The new redshift is outside that covered in a recent CO line search by Guil- 
loteau et at, (1999), highlighting the need for accurate redshifts for the obervation 
and interpretation of high-redshift line studies. 
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1 INTRODUCTION 

It is becoming increasingly clear that star-formation with 
its complementary photometric and spectroscopic signa- 
tures provides observers with a range of probes with which 
to study formation and evolution of galaxies in the early 
Universe. Particularly important to such studies are the 
millimetre-/submillimetre-/far-infrared- wavebands, a wave- 
length range in which clear fingerprints of both warm dust, 
heated by UV starlight from massive young stars, as well 



as atomic and molecular line emission, excited in regions of 
on-going s tarform ation, can be seen. Recent work by Wright 
& Reese ( 2000| ) has further illustrated the importance of 
this waveband, demonstrating that roughly half the radia- 
tion produced by starlight in galaxies is absorbed by dust 
and reradiated in the far- infrared/submillimetre- waveband. 
Likewise, COBE observations have shown that the cosmo- 
logical starformation rates inferred from UV- and optical 
studies (eg. Madau et ai, (1998)) are more than a factor 
of two lower than those required to produce the observed 
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Figure 1. A plot of the 850/xm flux vs. redshift for high-z objects with observed 850/xm fluxes and accurately determined rcdshifts. Stars 
denote detections of zjj4 radio-quiet quasars, with solid stars representing detections presented in this paper, and open stars sources 
taken from previous work. Squares denote detections of radio galaxies — unfilled taken from Archibald et al, 2000, and filled from the 
literature. Solid circles denote a selection of other high-z objects. Objects for which CO observations exist are ringed: solid rings - 
CO detections; dashed rings — CO non-detections. Lensing corrections are indicated by vertical dotted lines. The two curves trace the 
estimated flux that an Arp220-like object (based on a fit to observed FIR fluxes of Arp220) would have at different redshifts (solid line: 
Qm = 1; Ha. = 0.; dashed line: Qm = 0.3; Ca = 0-7) The dashed horizontal 3cr limit of this survey. 



far-infrared background, suggesting that there exists a pop- 
ulation of star-forming galaxies that are hidden at optical 
wavelengths. 

How can one use the (sub)millimetre-waveband to 



starformation potential. In addition spectral line measure- 
ments provide important kinematic and dynamic informa- 
tion. Whilst thermal emission from dust out to high red- 
shift can now be detected with relative ease due to the 



, „ ,. , ,, , ... ,. , new generation of sensitive multi-pixel millimetre (Kreysa 
tudy sj tartormation, both qualitatively and quantitatively, " i -n- — j 



the 



high redohift 
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conotrain tho root 



of the l (oub)millimotro wave omio 
frame spectral energy distribution - with more than 
one (sub)millimetre flux it is possible to evaluate a 
(sub) millimetre spectral index, thus confirming or refutin g 
a thermal o rigin o f the emission (eg. Hughes et al.( 1993), 
Isaak et al.( 1994)). With sufficient data, one can also de- 
termine a best-fit dust temperature, from which a bolo- 
metric (sub)millimetre-/far-infrared luminosity can be de- 
rived. Establishing the original energy source responsible 
for this reprocessed emiss ion is considerably more difficult 
(eg. Sanders et al, ( 1989)), as the continuum spectral sig- 



natures of dust-obscured AGN and massive starbursts are 
very similar, indeed, with current telescopes indistinguish- 
able at (sub)millimetre-wavelengths. If one assumes that 
some fraction of the (sub)millimetre/far-infrared luminosity 
arises from reprocessed UV-light from massive stars, then 
it is possible to infer starformation rates also. The pres- 
ence and detection of the molecular gas tracer CO, found 
ubiquitously in nearby starburst galaxies, points to future 



et al. 199£) and submillimetre-wave arrays (Holland et al. 



199£), detecting CO is more difficult, requiring a very ac- 
curate measure of the source redshift because of the limited 
bandwidth (Av ~ 1GHz), and thus small fractional band- 
width at 100-200 GHz (< 1%), of the current generation of 
millimetre-wave receivers. 



2 SURVEY DESIGN 
2.1 Choice of targets 

Recent attempts to study starformation in the high red- 
shift Universe have focused on identifying target objects in 
two quite distinct ways. The first of these - deep, unbiased 
surveys of blank sky - only recently became possible with 
the introduction of sensitive, multi-pixel (sub)millimetre- 
wave bolometer arrays. Surveys undertaken using SCUBA 



(Sub millimetre Continuum Bolometer Array (Holland et al. 



199£) ) at th e James Cler k Max well Telescope (eg. Hughes 
et al.( 1996| ), Barger et a/. (| 199fj| ), Smail et aL( |l997| ), Eales 
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et al.{ 1999)) have been particularly successful, identify- 



ing some tens of submillimctre-selected sources with 850/ira 
fluxes Ss50fj,m > 5m Jy. Classifying these sources, in partic- 
ular assigning redshifts, is complicated by (i) the limited 
accuracy with which one can determine the source positions 
using a single-dish telescope (ii) the large number of poten- 
tial, though very faint, optical counterparts revealed in deep 
observations of the submillimetre-wave source fields. As a re- 
sult it has been possible to begin to derive the starformation 
properties of only hand f ul of s ubmillimetre-selected sources 
(eg. Frayer et al.( 



1998 



2ooq )) 



In contrast, objects with known redshift provide a very 
convenient starting point from which to base searches for 
and studies of star-formation at high-redshift. The first 
submillimetre-wave observations of high redshift objects 



Barvai- 



were m ade in the early 1990s (Iciements et al. 199S 

An^onucci and Coleman 1992). Following on from the 



successful detection b oth of the (sub)millimetre continuum 



(by Clements et al, (| 1992|) ) and CO emission (by Brown 



& VandenBout (| 199lF and subsequently Solomon, Rad- 



ford and Downes ( 1992)) from the ultra-luminous infrared 
galaxy IRAS F10214+4724 at z ~ 2.3, (sub)millimetre stud- 
ies of the continuum and CO line emission at high redshift of 
a small sample of radio-quiet quasars taken from the Cam- 
bridge APM BR/BRI z ~ 4 QSO survey (jlrwin, McMa, 



Hazard_^99lh |5torrie-Lombardi et al. 1996 



Storric- 



hon 

Lombaidi et al. 2001)) , were made. At the time, these quasars 



were amongst the most distant known. With the c ombina- 
tion of the 800/^m/1100^ m measured with UKT14 (Duncan 
et al. 19^) on the JCMT ( |lsaak et al. 1994|) and the lSOO fim 
fluxes measured at IRAM-30m (iMcMahon et al. 19941), it 



was possible to show that the (sub)millimctre continuum 
spectral indices were consistent with thermal emission. The 
(sub)millimetre emission thus indicates the presence of mas- 
sive quantities of warm dust (Td ~ 50K) at redshifts of 
z ~ 4, synthesized, and possibly heated, by bursts of star- 
formation in the host galaxies. Large quantities of molecu- 
lar gas, Mh 2 > 10 Mq as traced by CO (eg . Ohta et al. ( 
199f^)7J)mont et al.( \ 1996a| ), Guilloteau et al. ( | 199S| )), were 
also detected in the host galaxies of a handful of quasars at 
z > 4. The detected gas masses are comparable to the total 
dynamical mass of the galaxy, a situation seen in galaxies 
undergoing their f irst bouts of star-formation (eg. Downes 
& Solomon (| 1998|)). 



2.2 Survey definition 

The number of objects from which both dust and CO emis- 
sion have been detected is small, and so it is currently not 
possible to differentiate between observed/derived proper- 
ties that are representative of objects at high redshift, and 
those that are extreme. An efficient way by which to in- 
crease the number of massive star-forming galaxies that can 
be readily studied needs to be identified. Shown in Figure 
^] is a plot of the observed 850/im flux vs. redshift for a se- 
lection of high-redshift objects with accurately determined 
redshifts. Looking at z J> 4: each source from which CO 
emission has been detected at greater than 5<j (see Guil- 
loteau et al., 1999 and references therein) has Ssso/jm > 
lOmJy, and conversely CO emission has been detected from 




0.2 



(225GHz) 



Figure 2. Plotted are the fraction of nights for the first 9 months 
of 2000 for which t 225GHz < r 225GHz' Note that the abscissa ex- 
tends to 0.2 only (and hence the maximum fraction of less than 1) 
as this represents the maximum t 225GHz f° r which 850fim obser- 
vations of the brightest sources can sensibly be made at JCMT (a 
zenith transmission > 45% ). Tau measurements were made us- 
ing the 225 GHz CSO tipping radiometer. Data have been taken 
from the CSO web-page, kindly provided by Ruisheng Peng. The 
dashed lines delineate the range of tau values during which the 
sky transmission was between ~ 60% — 70%, a T225GHZ range of 
0.08-0.12 (referred to as "grade 3" conditions at the JCMT). Such 
conditions occur for approximately 25% of the available night- 
time observing time. Note, the term "grade 1" is used to refer to 
occasions when the 850^tm zenith sky transmission is better than 
about 80% (T225 < 0.05), whilst "grade 2" conditions are those 
for which the zenith transmission is between about 70% and 80%. 

4 of the 5 objects with Sssom™ > lOmJy This suggests an 
observationally imposed cutoff for successful CO detections 
of Ss50fj,m ~ lOmJy. If we assume that by z ~ 5 (0.9(1.2) Gyr 
after recombination, using an Einstein De Sitter cosmology, 
with H = 50 kms -1 Mpc _1 and a A cosmology, Qm = 0.3, 
SIa = 0.7, H — 65 kms -1 Mpc _1 respectively) at least one 
massive burst of star-formation has already taken place, pro- 
ducing dust, then it is clear that it is more efficient with the 
current generation of radio telescopes to search first for the 
continuum dust emission from an object, following up secure 
detections with searches for CO. 

A 3cr detection of a source with Sssopm = lOmJy takes 
approximately 25 minutes (excluding pointing and calibra- 
tion overheads) using SCUBA/ JCMT under observing con- 
ditions for which the sky transparency is c. 65%. Such condi- 
tions occur ~ 25% of the time (Figure ^) and thus a survey 
for submillimetre-loud object at high redshift can slot well 
into periods during which the sky transmission is not high 
enough to undertake the high-priority projects requiring the 
highest possible sky transmission. Indeed, under the very 
best observing conditions, when the zenith transparency at 



* CO has not detected in BR B1117-11329, which has 
S^SSO^tm) = 13 ± 1 mjy, Buffey et al., in prep.) 
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Figure 3. Hubble diagram: absolute magnitude vs. redshift, of all known zj>4 radio-quiet quasars and a selection of z < 3 surveys. 
Absolute magnitudes have been calculated using a cosmology of H a = 50 kms -1 Mpc - , q a = 0.5, and magnitudes taken from the 
literature. The sources delineated by the solid box are those that define the optically-bright, z<;4 parent sample of the survey described 
in this paper, whilst the sources delineated bv the dashed box are those of the follow-up z ~ 2 sample. Black-hole masses have been 
calculated using the relationship given in Eq. fel and a bolomctric correction Lb i/Lg = 12 (Elvis et al. 1994) 



850^im can be as high as 90% , the required sensitivity can 
be reached in an elapsed time of around 12 minutes. Such 
observations become an inefficient use of telescope time since 
they become dominated by slew time and pointing checks. 
In contrast, the detection of associated molecular gas for a 
lOmJy source would currently take some tens of hours using 
the PdB Interferometer and the Owens Valley Millimeter 
Interferometer [j]. 

Samples of statistically significant size over a wide range 
of redshifts are crucial not only to studies of starformation 
and AGN activity at different epochs, but also to the deter- 
mination of the role of starformation activity and black-hole 



accretk n in the evolution of the Universe. 



In this paper we report on the first results of an 850/im 
program to search for massive star-bursts in the host galax- 



t If we assume that the submillimetre flux is thermal in ori- 
gin, then with assumptions about the physical properties of the 
dust: Td = 40 K, a dust spectral index of j3 = 2 (see Priddey 
& McMahon ( 2001) for a "fiducial" 2^4 quasar), a power- 
law dust emission given by Eq. |i| we derive a dust-mass esti- 
mate of around 10 9 M0for Ssso^m = lOmJy. Assuming a "typi- 
cal" gas-to-dust ratio of 500, and Mjj 2 /CO conversion factor of 
a = 4, the equivalent CO gas mass that one might expect is 
~few xlO 11 MQdetectable at millimetre wavelengths using either 
the Plateau de Bure interferometer or the Owens Valley Radio 
Observatory Interferometer in some 10s of hours 



ies of the most optically bright radio-quiet quasars at z <; 4 
using SCUBA at the JCMT. The aim of the survey is to 
identify a sample of ~10 submillimetre-bright sources, from 
which it should in principle be possible to detect molecu- 
lar gas emission using the current generation of millimetre 
interferometers. Based on extrapolations to 850^tm of ob- 
servations of radio-quiet quasars at z Jl 4 by Omont et al. 
( 1996b), we might expect to detect ~ 0.3 submillimetre- 
bright sources per hour (S^som™ ?=S lOmJy) under modest 
submillimetre-observing conditions (sky transparency of c. 
65% at 850/^m). We note that this detection rate is based on 
the sim plifying assumption that the sample of Omont et al. ( 



1996b) has the same distribution of source properties as any 



sample that we might define for the purposes of these obser- 
vations. Such newly identified submillimetre-bright sources, 
together with the 4 z ~ 4 well-studied objects from which 
CO has already been detected, and those for which upper 
limits to a gas mass have been obtained, will provide a sam- 
ple of statistically significant size (N > 3%/iV), with which it 
will be possible to begin to draw conclusions about massive 
star-formation and its role in the high redshift universe. 

We present here the interim results of our survey for a 
subset of 38 sources and discuss the implications of these 
current and future observations on studies of star-formation 



at high-redshift. In a companion paper, Omont et al. ( 2001) 
present 1.2mm observations of the same parent sample of 
zj;4 quasars. In a future paper(Priddey et al, in prep.) we 
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will present an analysis of the combined sample of more 
than 100 z<;4 quasars, including observations by others of 
different samples , which have now been observed at either 
or 1.2mm, 

We have adopted an Einstein de Sitter cosmology of 
H = 50 kms -1 Mpc" 1 , Q, — and, where useful, include 
a comparison with parameters derived using a A-cosmology 
with H = 65 kms" 1 Mpc" 1 , Q M = 0.3; fi A = 0.7. 



3 SAMPLE DEFINITION 
3.1 Optical properties 

The parent observational sample has been selected to in- 
clude the most luminous (rest-frame UV) 2^4 radio-quiet 
quasars known as of January, 2000. The quasar sample was 
primarily drawn from the APM multicolour photographic 
surveys based on APM scans of UKST plates by Irwin, 
McMahon & Hazard(1991), Storrie-Lombardi et a/.(1996, 
2001) and the Caltech DPOSS based on STSCI scans of 
POSS-II plates by Kennifick et aJ.(1995). The first step in 
the sample definition was to calculate the absolute mag- 
nitude in a consistent manner for all (~200) known zj>4 
quasars in order to select a sample on the basis of their 
rest-frame absolute B band magnitudes (see Figure [|). 

Absolute B band magnitudes {Mb) were evaluated us- 
ing observed R band magnitudes for all the quasars, assum- 
ing a spectral index a — —0.5 and using the relationship 



Me 



R - dm(z) + 2Mog(l + z) + k R (z) - 0.5, 



(1) 



where R is the observed i?-band(A ~ 6500A) mag- 
nitude and dm(z) is the distance modulus, evaluated 
at the source redshift. R band magnitudes were de- 
termined from APM scans of UKST OR (5900-6900A) 
or of POSS-1 E (6200-6800A) plates from APMCAT 
(www.ast.cam.ac.uk/~apmcat; McMahon & Irwin(1992)). 
Errors on these magnitudes are estimated to be or ~ 0.3. 
At z ~ 4, the R band starts to sample the rest-frame spec- 
trum near the prominent Lyman-a emission line; and at 
higher redshifts, the band suffers from the strong intergalac- 
tic absorption shortwards of Lyman-a. To compensate for 
both of these effects, a redshift- and filter-dependent cor- 
rection factor (Kr(z)) is invoked. Initially, for the purposes 
of sample selection, we adopted the empirical correction of 
Kennefick, Djorgovski & Meylan (1996). However, for the 
subsequent analysis, we derived our o wn correction, base d 
on the mean LBQS quasar spectrum (Francis et al. 1991), 
combined with the model of th e hydr ogen opacity of the in- 
tergalactic medium of Madau ( 1995| ), to obtain an average 
correction for each filter as a function of redshift. The cor- 
rection was compared with empirical values derived from all 
2 > 4 quasars for which we possess accurately flux-calibrated 
spectra (see Storrie-Lombardi et al., 1996, and Peroux et al., 
2001), and is found to follow the trend of the data well. The 
scatter superimposed on this trend is large, however, leading 
to an overall uncertainty in the magnitude aja ~ 0.5. This 
scatter reflects not only the idiosyncratic variations in the 
effects of emission and absorption features, but also an in- 
trinsic dispersion in the slope of the quasar continuum. For 
some of the quasars in the sample, a more direct measure 
of the continuum luminosity is available either through op- 
tical spectroscopy or through photometry in the K band, 



which, at z = 4, samples the rest-frame B band directly. 
However, for consistency and for the purposes of compari- 
son, we use the uniformly available R magnitudes to derive a 
homogeneous measure of luminosity. Issues pertaining to the 
derivation of accurate source magnitudes will be addressed 
more completely in a following paper. 

Objects already observed to a sensitivity equivalent 
or better than that of the survey as parts of other 
(sub)millimetre programs were excluded from the target 
list (eg. the 5 quasars observed by McMahon etal, 1999; 
BR B0019-1522, PSS J0134+3307, PSS J0747+4434, BR 
B1600+0729, BR B2237-0607), as were sources with tran- 
sit elevations of less than 60 degrees. An exception to this 
was PSS J1048+4407, previously observed and detected at 
1.25mm at the IRAM 30-m (Maoli et al., in prep), but not 
confirmed with 1.35mm observations a t the Plateau de Bure 
interferometer flGuilloteau et al. 1999| ). 



3.2 Radio properties 

The final source selection criterion was that of radio flux. 
We have used the NRAO VLA Survey Survey(NVSS) to 
measure or place upper limits on the radio emission from all 
z<;4 qu asars that were pote ntially in our sample. The NVSS 
survey (Condon et al. 1998) is a sensitive radio survey that 
was carried out with the VLA at 20cm in B-array, with 
a 5<t sensitivity of ~ 2. 5m Jy /beam. To determine which 
of the optically selected quasars are radio-loud, we cross- 
correlated the optical positions of our parent sample with 
sources in the NVSS catalogue. To assess the probability 
of a radio source being matched to one of the optically se- 
lected quasars by chance, the background source density was 
determined by counting the density of radio sources in the 
annulus with radius 100" to 1000", resulting in a value of 50 
deg" 2 . The rms positional error of the NVSS catalogue is 
given as 7"and thus, if one assumes that any radio source 
within 2l"is a potential match, one expects 0.005 chance 
associations per source, or a probability of 1 in 200 of in- 
correctly attributing radio emission to an optically selected 
quasar. Even if this radio emission were not to be from the 
quasar itself, a fraction of the emission would still be picked 
up by JCMT/IRAM main and error beams. 

Listed in Table [l] are all optically-selected zj> 4 quasars 
for which we find an NVSS source within a 30"radius - a 
total of 17 objects - along with their derived parameters. A 
detailed discussion of the radio properties of the sample is 
beyond the scope of this paper and will be presented else- 
where. All sources with Siaghz > 2.5mJy were excluded 
from the final SBQS target list. Based on a flux extrapola- 
tion of S oc v~ a , and a canonical spectral index of a = 0.7, 
we therefore expect the maximum synchrotron contamina- 
tion at 850/xm to be less than O.lmJy. Thus, it is reasonable 
to assume that the emission from any source that is detected 
at 850/j.m is likely to be thermal in origin. 

Our final target list comprised 76 sources with Mb < 
-27.5, augmented by 22 sources with -27.0 < Mb < -27.5 
which were included to improve sky coverage: 15 quasars 
were taken from the APM B/R and B/R/I survey ( BR and 



BRI sources : (Irwin, McMahon fc Hazard 1991 



Lombardi et al. 1996; Storrie-Lombardi et al. 2001)), 70 from 



the Palomar Sky Survey (PSS sources: (Djorgovski 200C)), 2 



Storrie- 
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Table 1. Optically and X-ray selected z ^> 4 quasars with 1.4GHz associations in NVSS. (Assuming a cosmology 
Ho = 50 kms" 1 Mpc - , qo = 0.5, and a spectral index for both radio and optical a = —0.5) 



Quasar name 


RA (J2000) 


Dec (J2000) 


2 


7 □ 


log P(r) | 


SlA 


„ n 


M B 


1 T El 

log 1/1.4 [1 


log Pi. 4 [| 












(false) 


mjy 


mjy 








PSS J0121+0347 4 


01 21 26.1 


+03 47 07 


4.13 


1.7 


-4.5 


78.6 


1.65 


-27.5 


34.69 


2.97 


PSS J0211+1107 


02 11 20.0 


+11 07 17 


3.99 


0.5 


-5.5 


44.1 


0.93 


-27.7 


34.41 


2.61 


PSS J0439-0207 


04 39 23.2 


-02 07 02 


4.40 


0.5 


-5.5 


43.8 


0.92 


-27.5 


34.49 


2.77 


PSS J2256+3230 


22 56 10.4 


+32 30 19 


4.04 


2.41 


-4.1 


14.7 


0.31 


-26.0 


33.95 


2.83 


BRI B0151-0025 A 


01 53 39.6 


-00 11 05 


4.20 


7.3 


-3.2 


3.0 


0.07 


—27.5 


33.28 


1.56 


BR J0234-1806 


02 34 55.1 


-18 06 09 


4.30 


0.7 


-5.2 


35.8 


0.75 


-27.5 


34.38 


2.66 


BR J0324-2918 


03 24 44.3 


-29 18 21 


4.62 


0.6 


-5.3 


236.8 


5.00 


-28.8 


35.26 


3.02 


BR J0355-3811 


03 55 04.9 


-38 11 42 


4.58 


13.4 


-2.6 


2.4 


0.05 


-29.4 


33.26 


0.78 


BR J0523-3345 


05 25 06.2 


-33 43 06 


4.40 


1.0 


-4.9 


188.7 


3.96 


—28.1 


35.12 


3.16 


BR J1053-0016 


10 53 20.4 


-00 16 49 


4.29 


2.6 


-4.1 


9.7 


0.2 


-27.0 


33.81 


2.29 


BR J1305-1420 


13 05 25.2 


-14 20 42 


4.04 


0.96 


-4.9 


20.6 


0.43 


-28.0 


34.09 


2.17 


SDSSp J0131+0052 s 


01 31 08.2 


+00 52 48 


4.19 


3.76 


-3.8 


4.2 


0.09 


-25.8 


33.43 


2.39 


SDSSp J0210-0018 s 


02 10 43.2 


-00 18 18 


4.77 


3.24 


-3.9 


11.5 


0.24 


-26.8 


33.97 


2.53 


SDSSp J0300+0032 s 


03 00 25.2 


+00 32 24 


4.19 


7.80 


-3.1 


7.1 


0.15 


-26.0 


33.70 


2.58 


SDSSp J1235-0003 c 


12 35 03.0 


-00 03 31.8 


4.69 


0.48 


-5.5 


19.5 


0.41 


-26.3 


34.19 


2.95 


PC J0027+0525 A 


00 29 49.97 


+05 42 04.4 


4.10 


8.47 


-3.0 


4.8 


0.10 


-24.8 


33.47 


2.83 


PC J2331+0216 A 


23 34 32.0 


+02 33 22 


4.09 


4.92 


-3.5 


2.7 


0.06 


-27.3 


33.22 


1.58 



For completeness all quasars in our parent sample for which potential radio counte rparts were found are li sted in 
the table. Those marked with A and B have already been found by Stern et of., J 200o|) and Fa n et al, ( 200l| ) 
to have NVSS counterparts, whilst the source lab elled c has a radio-counterpart (Fan et al. 2000) as found in the 
FIRST survey ( Becker, White and Helfand, 1995] ) 



a Separation between optical and nearest radio source positions 

b Logarithm of the probability that the optically selected quasar has been falsely associated with a radio source 
from the NVSS catalogue, separated by distance r 

c Extrapolated contribution by synchrotron emission to any observed 850fim flux, based on a spectral index, a, of 

0.7, S v tx v~ a 

d [L1.4] = ergs" 1 !^- 1 

e Ratio of radio (1.4 GHz) to optical (rest-frame B-band) luminosities 



from t he Palomar Survey (PC so urces : ( Schneider, Schmidt 
& Gun ^ 199l[ ISchncidcr 199l[ )), 1 from the VLA (VLA 
source: (3tern et al. 2000)) and 9 fro m the Sloan Digital 
Sky Survey (SDSSp sources: Fan et al.\ 199E) 



2000)) 



4 OBSERVATIONS AND DATA REDUCTION 
4.1 Observations 

Observations were made using the wide-band 850/xm fil- 
ter on SCUBA, during periods when the zenith sky trans- 
parency was ~ 60 — 75% (rsso,™ ~ 0.08 — rgso^m ~ 0.12, 
referred to as grade 3 weather in JCMT parlance). Sources 
were observed by the "displaced observer" with priority set 
by source availability and luminosity. The standard 850/im 
photometry observing mode was used - a 9-position jiggle 
superposed onto a 7 Hz/60"chop throw of the secondary 
mirror in the azimuthal direction, implemented to remove 
sky emission offsets. Gradients in sky emission were taken 
out by "nodding" the telescope - a physical movement of the 
telescope to the off-position every 18s. A typical observation 
consisted of sets of c. 30-50 18-second integrations, with the 
total number of integrations chosen to achieve a sensitivity of 
5*850,™ (1°") ~ 3.3mJy. Pointing was checked regularly, and 



found to be accurate to better than a few arc-seconds. As- 
tronomical seeing was monitored using the 5GHz CfA phase 
interferometer - the phase monitor was not in continual use 
however, when working, the data suggested that the seeing 
was not worse than 3"^J Sky transparency was monitored us- 
ing regular 850/ira sky-dips and data taken by the 225 GHz 
Caltech Submillimetre-Observatory (CSO) tau-meter. Flux 
calibration was achieved using the primary calibrators Mars 
and Uranus, as well as the secondary calibrators HLTau, 
IRC10216 |, OH231 and 16293 [[} Absolute fluxes for the 
planets were derived using the Starlink package, FLUXES. 
The average time to complete the observation of a single 
source, including pointing, calibration, sky-dip and slewing 
overheads, was over 1 hour. This is quite a bit longer than 



t We note that in the worst case, the combined pointing and 
seeing errors of 3" would result in a loss of point-source coupling 
of 40%, assuming a Gaussian beam profile with a FWHM of 14" 
§ IRC10216 was adopted as calibrator on one occasion only, due 
to its periodic variability. When used as a calibrator, the output 
flux of IRC10216 was at a maximum 

" See the JCMT web-page for the most recent calibrator fluxes 
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would be the case for a contiguous set of observations, re- 
flecting the start/stop nature of a back-up program. 

Here, we discuss observations of a subset of the target 
list: sources observed during the period February - July, 
2000. 



4.2 Data reduction 

Data we re reduced using the Sta rlink SURF reduction 
package (Jenness fc Lightfoot 200C), using the automated 



ORACDR pipeline (Jenness & Economou 2000) as well as a 



basic SURF-based pipeline written in-house prior to the re- 
lease of ORACDR. Two different schemes were implemented 
to determine an appropriate measure of the background sky- 
noise, the first using the mean of the inner ring of six pix- 
els surrounding the central on-source bolometer to define 
the temporal sky variations, and the second the median of 
all bolometers. In general, only small differences in the fi- 
nal fluxes and rms noise determinations were seen between 
the two methods for observations in which detections were 
made. There were, however, more significant differences in 
rms values determined for non-detections. In particular, the 
sky-noise correction derived from the inner-ring mean re- 
sulted in a non-zero offset in the distribution of the sig- 
nal/noise ratios for all off-source bolometers, indicating in- 
complete sky-removal. This was in contrast to the median- 
derived sky correction which we chose to adopt in the fi- 
nal analysis. All quasar data were clipped at the 20cr level 
prior to applying the sky-subtraction in order to eliminate 
transient spikes. Corrections for atmospheric extinction were 
made using sky opacities derived from 850^m sky-dips. In- 
dividual sets of integrations were gain-calibrated and then 
concatenated, with a 3cr despiking applied to remove outlier 
data points. 



5 RESULTS 

Source positions, redshifts and derived absolute B- 
magnitudes (Mb) are tabulated in Table 2, along with 
dates of observation, observed flux densities, rms values and 
signal-to-noise ratios. Note that calibration errors of ~20% 
have not been included in the error budget. 

A total of 38 quasars have been observed (see Table 2): 
8 sources were detected with Sgso/™ 10 m Jyi a t statisti- 
cal significances of better than 3a. 30 sources were observed 
but not detected: 19 sources with 3<r < 10 mjy, 8 sources 
with WmJy < 3<r < 12 mjy and 3 sources with 3<r > 12 
mjy (see Table 2). Achieving the pre-requisite survey sensi- 
tivity consistently proved to be difficult - the initial 3u ~ 10 
mjy sensitivity threshold was frequently both bettered (19 
non-detections) and not reached (If non-detections). 

The size of the observed sample is large enough not only 
for a significant number of submillimetre-bright sources to 
have been newly identified, but also to define a sample with 
which to start to assess the detection rate as well as de- 
tectability as a function of redshift and absolute magnitude. 
Here we adopt a 3<r ~12mJy cut-off (a limit which includes 
90% of sources observed in our program, Figure Q), noting 
that there are possibly as many as 3 sources which may have 



* 0.6 




0.4 - — - 



0.2 



2.5 3 3.5 

Observed RMS (mjy) 



Figure 4. A plot of the rms values obtained for all observa- 
tions reported in this paper. The histogram denotes the number 
of sources with a given rms, while the solid curve denotes the nor- 
malised cumulative number of sources with a given or lower rms. 
Note the dashed line delineates 90% of the objects observed, and 
defines the cutoff rms chosen for the statistical sample discussed 
in the text. 



Ss50/j,m ~ fOmJy, but for which the observations obtained 
were not deep enough to obtain a 3<r detection. Sources were 
initially included in our target list based on an approximate 
3cr ~ 10 mjy cutoff at 850/im - BR Bf 500+0824 and BR 
B 1302- 1404 were thus included in the JCMT sample in spite 
of having been observed at 1. 25mm at IRAM-30m down 
to an equivalent iSgso^m ~ 10 mjy. We include these two 
sources in our statistical sample which includes all sources 
for which the observed sensitivity is 3<r < 12 mjy - a total 
of 35 objects. 



6 DISCUSSION 

6.1 Newly detected submillimetre-bright sources 

To date, our survey has identified a total of 8 submillimetre- 
bright sources (where here we include PSS J1048+4407) 
with fluxes between ~ 10—25 mjy. This equates to a ~ 24% 
detection rate, which is com parable to the rate of ~ 30% 
seen by Omont et al. 



1996b) ,and based on 1.25 mm obser- 



vations of a sample of radio-quiet, optically selected APM 
quasars with z*;4. We stress here that the comparison be- 
tween the detection rates of the two different samples is 
valid ONLY if the underlying properties of the constituent 
sources are the same. We defer a more detailed analysis of 
this issue to another paper [J]. In each case, the measured 



II To make the comparison between the two different samples, 
observed at two different wavelengths, we use th e flux ratio of 
850Atm/1250/im derived by Priddey & McMahon ( 2001 ) for z^4 
quasars. At z ~ 4.2 this ratio is ~ 2.6, and thus we can compare 
the sample here with sources observed to la < 1.5 mjy at 1.25m 



© 2001 RAS, MNRAS 000, §-|l| 



K. G. Isaak, Robert S. Priddey, Richard G. McMahon et al. 



Table 2. Observational parameters of the observed z^>4 sample 
R RA Dec. R | M B ^ obs. 



(J2000) (J2000) 



date 



integ. | Flux ^ 
time (s) (mjy) 



rms 
(mjy) 



S/N 



PSS J0452+0355 4.38 04 52 51.5 +03 55 57 19.1 -27.6 7Feb00 2640 10.6 2.1 5.0 

PSS J0808+5215 4.45 08 08 49.4 +52 15 15 18.3 -28.7 4Feb00 2640 17.4 2.8 6.2 

PSS J1048+4407 1 4.38* 10 48 46.6 +44 07 13 19.6 -27.4 12Feb00/3Jul00 3940 12.0 2.2 5.3 

PSS J1057+4555 4.12 10 57 56.3 +45 55 53 16.5 -29.4 HFebOO 2640 19.2 2.8 6.9 

PSS J1248+3110 4.35 12 48 20.2 +31 10 44 18.9 -27.6 22Apr00 1340 12.7 3.4 3.8 

PSS J1418+4449 4.28 14 18 31.7 +44 49 38 17.6 -28.6 23Apr/3Jul00 4200 10.4 2.3 4.5 

PSS J1646+5514 4.04* 16 46 56.5 +55 14 46 17.1 -28.7 13/14.02/7Feb00 5240 9.5 2.5 3.8 

PSS J2322+1944 4.17 23 22 07.2 +19 44 23 17.7 -28.1 27Jul00 1600 22.5 2.5 8.9 



PSS J0007+2417 


4.05 


00 


07 


38.7 


+24 


17 


21 


18.3 


-27.6 


13Jun00 


1340 


1.2 


3.5 


0.3 


PSS J0014+3032 


4.47 


00 


14 


43.0 


+30 


32 


03 


18.7 


-28.1 


30May00 


1340 


10.6 


4.2 


2.5 


PSS J0133+0400 


4.15* 


01 


33 


40.3 


+04 


00 


59 


18.3 


-27.6 


13Jun00 


1340 


-1.0 


2.8 


-0.4 


BR B0300-0207 


4.25 


03 


02 


53.0 


-01 


56 


07 


18.5 


-27.6 


12Feb00 


3940 


-1.7 


2.8 


-0.6 


BR B0351-1034 


4.35* 


03 


53 


46.9 


-10 


25 


19 


18.6 


-27.8 


12Feb00 


2640 


-2.8 


3.0 


-0.9 


BR B0401-1711 


4.24 


04 


03 


56.6 


-17 


03 


24 


18.7 


-27.6 


HFebOO 


3940 


2.6 


3.4 


0.8 


PSS J0852+5045 


4.2 


08 


52 


27.3 


+50 


45 


11 


17.7 


-28.2 


9/12Feb00 


5240 


1.8 


2.2 


0.8 


PSS J0926+3055 


4.19 


09 


26 


36.3 


+30 


55 


05 


16.4 


-29.5 


4/llFebOO 


6540 


0.8 


1.9 


0.1 


BRI B0945-0411 


4.14 


09 


47 


49.6 


-04 


25 


15 


18.8 


-27.3 


5FebOO 


2640 


2.1 


2.2 


1.0 


PSS J0957+3308 


4.25 


09 


57 


44.5 


+33 


08 


20 


17.8 


-28.3 


HFeb/7MarOO 


5396 


3.7 


1.9 


2.0 


BRI B1013+0035 


4.40 


10 


15 


49.0 


+00 


20 


19 


18.8 


-27.9 


12Feb00 


2640 


5.5 


2.7 


2.0 


PSS J1026+3828 


4.18 


10 


26 


56.7 


+38 


28 


45 


18.1 


-27.8 


12Feb00 


2640 


5.2 


2.5 


2.1 


PSS J1058+1245 


4.33 


10 


58 


58.4 


+12 


45 


55 


17.6 


-28.8 


12Feb00 


2640 


0.00 


2.2 


0.0 


BRI B1110+0106 


3.92* 


11 


12 


46.3 


+00 


49 


58 


18.3 


-27.1 


5Apr00 


1340 


7.1 


3.4 


2.1 


SDSSp J1226+0059 


4.25 


12 


26 


00.7 


+00 


59 


21 


18.9 


-27.0 


14Feb00 


2640 


-0.06 


3.3 


0.0 


BR B1302-1404 


4.00* 


13 


05 


25.2 


-14 


20 


12 


18.6 


-27.4 


23Apr00 


1600 


1.0 


3.2 


0.3 


PSS J1315+2924 


4.18 


13 


15 


39.8 


+29 


24 


39 


18.5 


-27.4 


22Apr00 


1340 


4.1 


3.3 


1.3 


PSS J1326+0743 


4.17 


13 


26 


11.8 


+07 


43 


58 


17.5 


-28.4 


3Jul00 


1340 


-0.3 


3.9 


-0.1 


PSS J1347+4956 


4.46 


13 


47 


43.3 


+49 


56 


21 


18.5 


-28.3 


22Apr00 


1470 


8.5 


3.9 


2.2 


SDSSp J1413-0049 


4.14 


14 


13 


32.4 


-00 


19 


10 


19.1 


-27.0 


14Feb00 


2640 


-1.0 


3.2 


-0.3 


PSS J1432+3940 


4.28 


14 


.32 


24.8 


+39 


10 


21 


18.0 


-28.2 


3Jul00 


2640 


-1.4 


3.3 


-0.4 


BR B1500+0824 


3.94* 


15 


02 


45.4 


+08 


13 


06 


19.3 


-26.3 


14Feb00 


3940 


6.2 


3.0 


2.1 


PSS J2154+0335 


4.36 


21 


54 


06.7 


+03 


35 


39 


19.0 


-27.5 


5Apr00 


1860 


2.6 


2.6 


1.0 


PSS J2155+1358 


4.26* 


21 


55 


02.2 


+13 


58 


26 


18.0 


-28.1 


5Apr00 


1860 


0.8 


3.1 


0.1 


PSS J2203+1824 


4.38 


22 


03 


43.4 


+18 


28 


13 


18.0 


-28.7 


30May00 


1340 


3.0 


4.2 


0.7 


BR B2212-1626 


3.99* 


22 


15 


27.3 


-16 


11 


33 


18.1 


-27.9 


30May00 


1340 


7.9 


4.5 


1.7 


PSS J2238+2603 


4.03 


22 


38 


41.6 


+26 


03 


45 


17.0 


-28.9 


30May00 


1340 


2.2 


3.9 


0.6 


PSS J2241+1352 


4.44* 


22 


41 


47.8 


+13 


52 


02 


19.1 


-28.0 


5Apr/30MayOO 


2250 


-0.2 


3.0 


-0.05 


PSS J2323+2758 


4.18 


23 


23 


40.9 


+27 


58 


00 


18.6 


-27.3 


30May00 


1340 


-2.8 


4.0 


-0.7 


PSS J2344+0342 


4.24* 


23 


41 


03.2 


+03 


42 


26 


18.2 


-28.1 


30May00 


1340 


11.3 


4.0 


2.8 



a BR/BRI source posi tions taken from the APM catalogue; PSS source parameters have been taken from the 
Palomar Sky Survey ('([Diorgovski 2000[)) SDSSp source positions have been taken from the Sloan Digital Sky 



Survey (( [Fan et al. 1999 ; |Fan et al. 2000[ )) 

b Palomar Sky Survey redshifts have been taken from Djorgovski ( 200C), with the exception of sources marked *, 
where redshifts have been taken from Peroux et al. and Storrie-Lombardi et al. 

c R magnitudes determined from APM scans of the UKST OR (5900-6900A) or POSS-1 E (6200-6800A) plates. 
Errors are estimated to be <tr ~ 0.3 

d See text for an explanation of the method used to evaluate the absolute source magnitudes 

e Integration times have been calculated based on no: integrations x 18s/integration x 1.5, where both chopping 
and a 50% observational SCUBA overhead has been allowed for, This does not include any overhead for calibration 
or pointing. 

f Where a 3<r detection was not been achieved, the tabulated flux is given as the on-source signal only 
9 The position reported here is that taken from the Djorgovski web page. A check of the position using the APM 
catalogue after the JCMT SCUBA observations had been made revealed that the published position is out by 2". 
The correct position is: (J2000) 10 48 46.64 +44 07 10.9 

h The flux reported is the weighted sum of 4 sets of integrations that differed considerably in statistical significance, 
ranging from 3<r and better detections to a near-zero on-source flux. 
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fluxes are more than a factor of 100 in excess of what might 
be expected based on the NVSS radio limits. Thus, the sin- 
gle submillimetre flux measurements, combined with the a- 
priori radio limit of our parent sample, suggest strongly that 
we are detecting submillimetre emission from these sources 
that is thermal in origin. 

How do the fluxes reported here compare with those 
measured for other z<;4 radio-quiet quasars? The brightest 
radio-quiet quasar zjj;4 host galaxy observed both at sub- 
millimetre and mi llimetr e wavelengths re mains BR B1202- 
0725 (Isaak et al. ( | 1994j ), McMahon et al.( \ 1994| ), Buffey et 
al, in prep) - with an 850/im flux of almost greater than 
twice that of the brightest sources detected in this survey, 
PSS J0808+5215, PSS J1057+4555 and PSS J2322+1944. 
Arc-second resolution ( c. 1-2") imaging with the Pl ateau 
de Bure interferometer (PdB) by Omont et al.( 1996b ) have 
shown BR B1202-0725 to be a composite source, made up of 
two comparable gas-rich and possibly interacting galaxies. 
If we assume that the 850/^m flux distribution is similar to 
that of the 1.25mm flux, split equally between the two ob- 
jects, then BR B1202— 0725 is no longer such an exceptional 
object, and has a flux comparable to that of the brightest 
sources present here. With the further assumption that PSS 
J0808+5215, PSS J1057+4555 and PSS J2322+1944 are all 
single sources, there is a suggestion of a maximum observed 
flux of Sg50/!tm ~ 20 mjy. The origin of such a maximum 
is intriguing. A possibility is that the starformation is self- 
regulated, with the large number of supernovae events in 
the massive star-burst disrupting the local ISM and block- 
ing star-formation. We are currently obtaining both deep 
radio observations with the VLA and millimetre observa- 
tions with PdB, with which we will be able to establish the 
degree of extension of the sources. 

We note that the sources we detect here are as bright, 
if not brighter than those detected in blank-sky and cluster- 
lensing surveys. Interestingly, the two brightest sources de- 



tected in the cluster lensing surveys (Ivison et al. 1998 
Knudsen et al. 2000) have 850/im fluxes of just over 
3(850^) ~ 20m Jy, and are both AGN at z ~ 2. For com- 
parison, the measured fluxes of the brightest and median 
source in the Canada-UK-De ep-Sky-Survey are ~ 8mJy and 
~ 4.5mJy respectively (eg. flEales et al. 1999j )). The 200- 
square arcminute survey "8mJy survey" (3a <8mJy (Dun- 
lop 2000, Scott et al, in prep.) has yielded 24 sources, de- 
tected at statistically significant levels, with 10 brighter than 
lOmJy, and the brightest with a flux of 15mJy. We note here 
that others have also started similar surveys to o ur own at 
1.25mm using the IRAM-30m jCarilli et aLWQla ), but leave 
a comparison between surveys to a companion paper. 



6.2 PSS J1048+4407 

In addition to identifying new submillimetre-bright sources, 
we detected PSS J1048+4407, a source for which observa- 
tions at ~ 1.3 mm with IRAM-PdB and the IRAM-30m 
were in conflict. 



6.2.1 Single-dish vs. inter ferometric observations 

As reported by flGuilloteau et al. 1999| ), PSS J1048+4407 
was not detected at 1.35mm using the PdB interferome- 
ter (Si. 35 mm = 0.25 ± 0.68 mjy) in the CD configura- 
tion, but was tentatively detected by Maoli et al (in prep.) 
1.25mm at the IRAM-30m (Si.asmm ~ 3 ± lmjy). Our re- 
sult, S$50fim = 12.0±3mJy, is consistent, to within the error- 
bars, with the 30-m single-dish measurement if we assume 
a ratio of ffgso^m / Si.2Smm ~ 2.5, as derived by Priddey & 



McMahon ( 2001). Neither of the two single-dish measure- 
ments, however, are consistent with the interferometer ob- 
servations. Pointing errors in the interferometer map and the 
effects of anomalous refraction would need to be much larger 
than realistically possible to explain the non-detection with 
the interferometer, given that at ~ 220GHz the primary 
beam is ~ 27". It is more likely that the host galaxy of PSS 
J1048+4407 is spatially extended on the scales greater than 
the synthesized beam of the interferometer (~ 3"x 2"). 

To resolve out a substantial fraction of the millimetre 
flux, the spatial extent of the emission does not need to be 
much greater than the size of the beam, depending more crit- 
ically on the projected baselines/effective U-V coverage of 
the observations and the source extension relative to the po- 
sition angle of the beam. At a redshift of z ~ 4.4, l"equates 
to a linear size of c. 6kpc To resolve out emission over arc- 
second scales suggests that the star-burst activity, as traced 
by the millimetre emission, is taking place over the whole 
galaxy rat her than confined to nuclear regions as observe d 
in Arp220 ( |Downes fc Solomon 1998j |Sakomoto et al. 1999] ), 
Arpl93 and Mrk273 ( pownes fc Solomon 1998[ ). ThisTsin 
contrast to BRI B0952-0115 and BR B1335-0417 where a 
comparison between the single-dish and interferometric ob- 
servations suggests that in each case a significant fraction of 
the millimetre emission is unresolved, but similar to the case 
of BR1117-1329 and BR1144-0723, as discussed in Omont 



et al, ( 2001) 



6.2.2 Updated redshift 



Also reported in Guilloteau et al., ( 1999) is the non- 
detection of CO (5-4) line emission using the PdB interfer- 
ometer. This is interesting in light of the 850/im detection 
presented here, as the source is one of two from five with 
a previously inferred Sssoum > lOmJy for which a search 
for CO emission failed. A reassessment of the source red- 
shift based on a new o ptical spectrum taken at the WHT 
by Peroux et al. ( 2001) (Figure fh reveals that the CO ob- 



by Omont et al.( 19961). In evaluating our detection rate we 
have excluded BR B1500+0824 and BR B1302-1404 because of 
the common non-detection to the IRAM-30m sample 



servations were centred at an erroneous redshift. Using the 
combination of the Lya, CIV and the SilV+OIV emission 
lines, we derive a redshift of z — 4.38, in contrast to tha t 
derived from the Lya line only by Kennefick et al.( 1995), 
z = 4.45, a nd the redshift at whic h a search for CO emission 
was made ( Guilloteau et al. 1999] ) . The difference in redshift 
is significant - much greater than that covered by the inter- 
ferometer correlators. Our average redshift may, however, 
be biased high as it contains the measurement based on the 
Lyman-a line which is observed to be systematically higher 
than the other two emission lines(CIV, SilV+OIV). It is 
possible that the Lyman-a line profile is affected by absorp- 
tion due to the NV(1240A) line which would bias its peak 
redward. Using the CIV line alone, we derive a systemic red- 
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Figure 5. An optical spectrum of PSS J1048+4408, taken at a resolution of 5A using the ISIS spectrograph on the WHT by Peroux et 
al. . On the basis of this spectrum, we classify this quasar as a BAL, not possible with previous spectra because of the limited spectral 
coverage. The presence of broad absorption lines makes an accurate determination of the host galaxy redshift difficult, however we derive 
redshifts of: z = 4.422 (Lya{1216A), z = 4.367 (Si+O (1400A) and z = 4.354 (CIV(1539A)). From an average of these three lines we 
derive a source redshift of z = 4.38, in contrast to z = 4.45 as first derived by Kennefick et al. We note that even this redshift estimate 
could be biased high, as explained in the text 



shift of z — 4.360, correcting for the systemic shift from rest 
for the line as derived by Tytler & Fan (1992). 

Given the discrepancy between the observed and newly 
derived redshift we suggest that the failure to detect CO 
emission is due to incorrect redshift, rather than an anoma- 
lous gas-to-dust ratio. This serves as a note of caution - 
given the limited bandwidth of the current generation of 
millimetre and submillimetre telescopes, the non-detection 
of any spectral line emission should be considered both in 
an astronomical context and, as important, in terms of a 
possible error in redshift. 



6.3 Inferred astrophysical properties 

Shown in Table 2 are physical properties of the detected 
sources derived from the optical luminosities and submil- 
limetre fluxes reported here. We have assumed that the 
measured submillimetre flux arises from thermal emission 
from warm dust, and that the dust is heated solely by the 
UV flux from massive star- formation. The former can be 
justified, at least in part, by the radio-selection criterion 
that we used to define the parent sample. In addition, we 
will be able to combine the 850/im fluxes measured here 
with 1.25mm fluxes that hav e been measured in a paral- 
lel 1.25mm MAMBO survey ( |Omont et al. 200l| ) running 
at the IRAM-30m (Priddey et al., in prep.). The latter as- 



sumption is more difficult to justify, as it is possible that a 
fraction of the rest-frame UV absorbed by the dust comes 
from the central AGN. The size of the AGN contribution is 
difficult to estima te as it depends on the source geometry 
(eg. Sanders et al.( 1989)), about which we have little knowl- 



edge. Observations of local, ultra-luminous infrared galaxies 
(ULIRGs, with Lfir > 10 12 Lq), arguably the most likely 
low-redshift counterparts of the quasar host galaxies, have 
been interpreted as being suggestive that between 70-80% 
of ULIRGS are powered predominantly by star-burst activ- 
ity, with at least half the sources observed probably have 
simultaneous star-burst and AGN activity on the ph ysical 
size scales of a few kpcs (eg. Genzel et al, ( 1998)). An 
additional unknown at this point is the extent of gravita- 
tional lensing and thus source magnification. The incidence 
of gravitational lensing amongst high lum inosity quasars i s 
around 1% for samples of bright quasars (Kochanek 1993). 
Studies of previous samples have, however, a median redshift 
of only z ~ 2. At z ~ 4 the probability of a quasar being 
lensed by a foreground galaxy is a facto r of ~ 2 higher due 



to the larger optica l depth to lensing (Turner et al. 1984 
Fukugita et al. 1992). This qualitative argument is in rough 
agreement with the empirical result that only one of the 49 
2 > 4 APM BRI quasars, the z = 4.43 APM quasar BRI 
B0952-0115 with Seso^m = 14mJy (Buffey et al, in prep) 
is known to be gravitationally lensed. It therefore seems un- 
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Table 3. Derived properties of the detected (submillimetre-loud) zj>4 radio-quiet quasars sources - see Section 6.3 for the 
definitions used in deriving source properties 



Source name z M B uL B Mbh M acc M sph S$ 5 o M D L FIR SFR t* 

10 13 L Q 1O 9 M ^ IO^Mq mjy 10 9 M 10 13 L J^£ 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


PSS J0452+0355 


4.38 


-27.6 


0.9 


3.6 


80 


7.3 


10.6 


0.9 


1.2 


1200 


0.6 


PSS J0808+5215 


4.45 


-28.7 


2.5 


10.0 


220 


20 


17.4 


1.5 


2.0 


2000 


1.0 


PSS J1048+4407 


4.38 


-27.4 


0.8 


3.0 


70 


6.0 


12.0 


1.1 


1.3 


1300 


0.5 


PSS J1057+4555 


4.12 


-29.4 


2.8 


19.0 


120 


38 


19.2 


1.8 


2.3 


2300 


1.7 


PSS J1248+3110 


4.35 


-27.6 


0.9 


3.6 


80 


7.3 


12.7 


1.1 


1.4 


1400 


0.5 


PSS J1418+4449 


4.28 


-28.6 


2.3 


9.1 


200 


18 


10.4 


0.9 


1.2 


1200 


1.5 


PSS J1646+5514 


4.35 


-28.7 


2.5 


10.0 


220 


20 


9.45 


0.8 


1.1 


1100 


1.9 


PSS J2322+1944 


4.17 


-28.1 


1.4 


5.8 


130 


12 


22.5 


2.0 


2.6 


2600 


0.4 


median z^4 SBQS Fl 
Fl 


4.35 


-28.4 


1.8 


7.4 


165 


15 


12.4 


1.1 


1.4 


1350 


1.1 


4.35 


-28.7 


2.5 


10.2 


230 


21 


12.4 


1.5 


2.0 


1950 


1.1 



a The quantities and parameters evaluated for the median quasar are based on the median observables, from which param- 
eters are derived. 

b Median quasar properties evaluated as above, based on a A cosmology, Ca = 0.7, Qm = 0.3, and H = 65kms — 1 Mpc . 
For comparison the ratio of the square of the luminositiy distance under the two cosmologies at z = 4.35 is 1.38 



likely that a large number of the submillimetre-bright source 
detections reported here are lensed. High resolution optical 
imaging of the whole sample would be needed to confirm 
this. At present we have no evidence to support or refute 
the possibility of magnification, and thus assume that there 



6. 3. 1 Astrophysical properties inferred from the optical 
AGN luminosity 

An advantage of targeting quasars for a study of high- 
redshift (sub)millimetre sources is that we are able to use 
the optical luminosity of the AGN to constrain the mass of 
the black hole powering the AGN, from which we can esti- 
mate the mass of the surrounding host galaxy. If we assume 
a bolometric correction from the B band of L^ \/Lb = 12 



(Elvis et al. 1994), and that the black hole is accreting at the 
Eddington rate, then there exists a simple relation between 
black hole mass, Mb, and its bolometric luminosity, Lboi: 



4nGm p cM hh 13 Mbh t 

^bol = = O X 11) -t n(\ K T ^©J 



10 9 M e 



(2) 



where o~t is the Thomson scattering cross-section and m p 
the proton mass. 

The accretion rate can be written in terms of the lumi- 
nosity and an efficiency par ameter e, fo r which a conservative 
estimate is eo.i ((Rees 1984)): 



Mace = = 



22 M B 



ec 2 e .il0 9 M c 



■M 



©yr 



(3) 



While Equation g holds, the growth of black hole mass 
and luminosity is exponential, with an e-folding time-scale 
face = Mbh/ Mlcc = ep.i x 4.5 x 10 7 yr. Studies in the lo - 
cal Universe (Vlagorrian et al. 1998; Gebhardt et al. 200C) 



have shown that there exists a linear relation between the 
mass of non-active black holes and the stellar mass of the 
surrounding galactic bulge at the current epoch. If we as- 
sume that there is some imprint of this relationship at early 



epochs (however, see Omont et al. ( 2001) and Priddey et 
al. (in prep)), then we can obtain a measure of the mass of 
the final (ie. 2 = 0) quasar host galaxy (Column 7 of Ta- 
ble 3) using a value of Mbh ~ 0.005M sp h from Gebhardt 
et al. ( 200C). It is clearly possible that at z = 4, the ra- 



tio of the black hole to spheroid mass is different from the 
local value. For example Kauffmann & Haehnelt (2000) pre- 
dict that the host masses will less massive by a factor of 
between 5-10 at z = 2 as compared to z = 0. Work by Rix 
et al (1999) on the magnitudes of lensed host galaxies of 
lensed quasars also suggests that the masses of quasar hosts 
are low at high redshift however although this is based on 
the observed rest-frame B band magnitudes, these observa- 
tions could be affected by dust obscuration, and hence the 
inferred masses may be underestimated. In the long term, 
one would hope that it will be possible to determine the dy- 
namical masses of quasar hosts at high redshifts so that the 
issue can be studied more directly. 

These estimates represent lower limits to the mass of the 
central engine, as we observe the quasars while they are still 
accreting. It seems unlikely for a number of reasons, how- 
ever, that the masses grow much larger than this. Firstly, the 
inferred masses are large, comparable to the most massive 
non-active blac k holes (;$ 10 9,5 Mq) a t the cores of nearby 
galactic bulges (Gebhardt et al. 2000). Secondly, the accre- 



tion rates required to sustain these luminosities (Table 3, 
Column 6) are extreme — and encroaching upon (~ 10% of) 
the star-formation rates inferred from the far-infrared lumi- 
nosity (Column 11). Thus, the gas supplies in the galaxy 
are being rapidly diminished, at a rate too high to sus- 
tain either processes for much longer than 1 Gyr - indeed 
Priddey & McMahon determine a median gas co nsumption 
timescale of 0.1 Gyr ([Priddey fc McMahon 2001]). Although 



we have no knowledge of the precise accretion history of the 
AGN, we can derive a time taken to assemble the inferred 

black masses of around 0.5 1 Gyr if we assume that the 

accretion rate is Eddington-limited throughout the growth 
of the black hole and that the seeds from which the black 
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holes form are M 
1998). 



10 3 " 6 M (Haehnelt, Natarajan & Rees, 



6.3.2 Astrophysical properties inferred from the 
submillimetre fluxes 

Turning to the host galaxy properties derived from submil- 
limetre flux observed at frequency v — f res t/{l + z), the 
dust mass given in Column 9 of Table 3 is determined using 



M D 



S{y )D\ 



(4) 



(l+z)Kz>(v r )B u (v T ,T D )' 

where we assume that the emissivity k is a power-law func- 
tion of frequency, and we adopt the normalisation of Hilde 
brand (1983): 



l2.4THz/ 



(•») 



,2.4THz, 

where ko = 18.75 cm 2 g _i is the value determined at A = 
125/im. For the dust temperature and emissivity index, we 
use T D = 4 0K and = 2.0, derived by Priddey & McMa- 



hon ( 2001) from a composite submillimetre spectrum of 
z^;4 quasars. The far-infrared luminosity is calculated by 
integrating under the thermal, grey-body spectral energy 
distribution, 



Lfir = 2nhc 



4+/3 



r(4 + /?)C(4 + /3), 



(6) 



V he J 

where again we assume T = 40K, j3 = 2.0, V is the Gamma- 
function and £ is the Riemann-Zeta function. If we assume 
that all the energy re-radiated by the dust originates from 
stars, then the far-infrared luminosity can be used to esti- 
mate the current star-formation rate: 
£fir 



M* = * x 



10 10 L c 



M 



©yr 



(7) 



The exact value of \E' depends sensitively upon the star- 
formation history of the host galaxy, the stellar Initial Mass 
Function (IMF) and the fraction of starlight absorbed by the 
dust cloud. We adopt the simplifying assumptions that the 
starformation rate is constant and that the cloud-covering 
factor is unity. Then, for a Salpeter IMF with upper and 
lower mass cut-offs of 100 M and 0.3 MQrespectively, ^ sa 1 
and is approximately constant for ages <; 10 8 yr. The to- 
tal time required to form the final stellar mass (Column 
7) inferred from the black hole mass is given in Column 
12 of Table 3, and is roughly 1 Gyr. The starformation his- 
tory derived from the chemical/dynamical model of ellipti - 



cal galaxy evolution developed by Friaga & Terlevich ( 199£ ) 
has just such a form — a constant star formation rate lasting 
for around 1 Gyr. Observational support for long periods of 
sustai ned hig h starformation rates is reported in Abraham 
et al. ( 1999) (see Figure 6). F or an alternative interpreta- 
tion, however, see Omont et al. ( 2001) who argue for shorter 
burst timescales. 

We have two time-scales describing the fate of gas in the 
evolving host galaxy: a time-scale for the conversion of gas 
into stars, ~ 1 Gyr, or gas accretion onto a super-massive 
central black hole with an accretion time of the same order. 
This may be suggestive of a possible synchrony between for- 
mation of the AGN and its host galaxy, however we leave 
a more detailed and speculative discussion of the inferred 
quasar and host galaxy properties to a future paper. 



Observed sources 
Delected sources 
Undetected sourc 




Absolute magnitude (M B ) 
Figure 6. Cumulative distribution of the fraction of detections 
(heavy solid line), non-detections (dashed line) and observations 
(light solid line) as a function of absolute magnitude. There is 
some indication of a difference in the cumulative distribution, 
however using the K-S test, it can be shown that the null hy- 
pothesis that the magnitudes of the detected sources and the 
non-detected sources are similarly distributed cannot be rejected 
at the 90% level. 



Throughout, we have assumed that the dominant heat- 
ing source is starformation in the quasar host galaxy. At 
present, however, it is not possible to determine conclu- 
sively whether this is indeed the case. If instead, we were 
to that assume the AGN contribution to the dust heating 
cannot be neglected, then (a) for a fixed dust temperature 
the inferred starformation rate would fall, with a correspond- 
ing increase in the formation time for the massive spheroid 
(b) a hotter dust component heated directly by the AGN 
would contribute significantly to the overall far-infrared lu- 
minosity, and thus inaddition, smaller dust masses as well as 
lower starformation rates and longer timescales would be in- 
ferred. Observations of samples, such as this SBQS sample, 
at shorter wavelengths using SIRTF will reveal any warmer 
dust component. 



6.4 Optical/submillimetre correlations 

With a sample of 38 observed sources, it becomes mean- 
ingful to start to make both qualitative and quantitative 
assessments of the trends in the detection rates with source 
redshift and magnitude. Shown in Figures ^ and |t] are the 
cumulative magnitude and redshift distributions for the ob- 
served, detected and non-detected subsamples. To assess 
quantitatively the degree of similarity between the distribu- 
t ions w e use the Kolmogorov- Smirnoff (K-S) test (eg. Barlow 
( | 1998; )). In Figure | we note that there is some indication of 
a difference between the detected and non-detected cumu- 
lative magnitude distributions however, formally, the differ- 
ence is not statistically significant, and so we cannot reject 
the null hypothesis that the two samples are drawn from the 
same population. The magnitude and redshift distribution 
of the two samples are shown in a slightly different way in 
Figures ^ and |y. 

Plotted in Figure [u] are the 850/im fluxes and absolute 
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Observed sources 
Detected sources 
Undetected sources 



Redshift 

Figure 7. Cumulative distribution of the fraction of detections 
(heavy solid line), non-detections (dashed line) and observations 
(light solid line) as a function of optical redshift. Again, using the 
K-S test, we cannot reject the null hypothesis that the redshift 
distributions of the detected and non-detected sources are the 



magnitudes for the detections presented in this paper, as 
well as upper limits (1.65a (90% confidence level) + "on- 
source signal") for non-detections. Error bars for both de- 
tections and non detections are la. It is difficult to make 
a quantitative assessment of the significance of the opti- 
cal/submillimetre magnitude correlation, given the signifi- 
cant fraction of upper limits in the sample. There are two 
competing physical processes that would produce a statisti- 
cally significant correlation between detected submillimetre 
flux and absolute quasar luminosity. Firstly, the most lu- 
minous quasars reside in the most massive host galaxies, in 
which there is likely to be ongoing star-formation - the larger 
the host galaxy, the greater the potential for massive bursts 
of star-formation, and thus more dust and more UV photon 
flux with which to heat the dust. More quantitatively, using 
Equations ^, ^ and 0, 



0.1 M, 
*eo.i Mice 



(8) 



If we observe these objects just before the end of their accre- 
tion lifetime, when the accretion rate encroache s upon the 
star for mation r ate, M acc ~ M* =>■ Lfir ~ vLb (Priddey & 



McMaf on 2001). Secondly, the more luminous the central 



quasar, the greater its UV-output, and thus the larger the 
power that may be absorbed by dust, and re-emitted in the 
rest-frame far-infrared. 

A plot of the source magnitude vs. redshift for detec- 
tions and non-detections also fails to reveal a preferential 
range of redshifts/magnitudes for which the number of de- 
tections is high. It must be cautioned, however, that the 
number of sources in any one bin is small, so the overall 
statistical significance of this result is low. 

In addition, a potential systematic bias arises from the 
fact that the absolute B-band magnitudes have been de- 
rived from the measured i?-band APM magnitudes. As dis- 
cussed in Section 3.1, calculating an optical luminosity from 
an R magnitude is complicated for sources at z <; 4 because 



Absolute Magnitude (M D ) 
Figure 8. Histograms of the absolute magnitudes (Mb) of the 
observed sample (dotted shading) and the detected sample (solid 
horizontal shading). 



Figure 9. Histograms of the redshift distribution of the observed 
(dotted shading) and the detected sample (solid horizontal shad- 
ing). 



the band is contaminated by redshifted Lyman-a emission, 
and at higher redshifts, intergalactic absorption. We intro- 
duced a correction to account these effects, and found that 
it worked well on average, but that the scatter about the 
average was <tm ~ 0.5. Obtaining optical spectra and/or K- 
band photometry for all the targets may ultimately help to 
alleviate this. 



6.5 Non-detections — an assessment of the 

contribution of the z ~ 4 population to the 
submillimetre-wave background 

Shown in Figure [H] is the signal-to-noise distribution of the 
non-detections (solid line) - this has a clear non-zero mean. 
More quantitatively, the weighted mean of fluxes of the ob- 
served non-detections is 2.0±0.6mJy, a quantity that we can 
evaluate because our survey has a sensitivity that is consis- 
tent between sources to within a factor of ~ 2. The mean 
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-27 -28 -29 -30 

Absolute Magnitude (M B ) 
Figure 10. A plot of 850/xm flux vs. Mb- Filled circles denote 
source detections reported in this paper, while crosses denote up- 
per limits to non-detections (where plotted flux = "signal" + 
1.645<t (90% confidence)). The flux error bars denote la errors, 
as given in Table 2, whilst magnitude error bars, ctmb ±0.5, are 
a measure of the combined systematic and random error in the 
measurement and derivation of the absolute magnitudes. 



o 
c 
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S/N ratio 

Figure 11. Histogram of the S/N ratios of all non-detections 
in the current sample: heavy-lined histogram denotes S/N ratios 
evaluated for the on-source pixel. Superposed (light-line) is the 
outline of the S/N distribution for all functional off-source pixels, 
over which is plotted a Gaussian approximation to the off-source 
distribution. Note that the mean of the off-source distribution is 
zero, and the distribution itself very nearly Gaussian. This sug- 
gests that the sky-offsets are being removed correctly. 



can therefore be used to provide a statistical measure of 
the submillimetre flux of the underlying radio-quiet quasar 
host galaxy population. The derived value is consistent with 
the earlier estimate of ~ 3mJy based on the dee per obser- 
vations of a sample 6 zjj4 radio-quiet quasars ( [McMahon 
et al. 1999), and is also around the value that one might 



and a mean flux as determined above, then the contribution 
of the quasar host galaxies to the 850/^m background is neg- 
ligible. 



7 FOLLOWUP WORK 

Whilst useful for identifying the most submillimetre-bright 
sources, the detection of objects at a single wavelength 
clearly provides only limited insight into the detailed phys- 
ical properties of the individual sources. Given the strong 
dependence of the inferred far-infrared luminosity on dust 
temperature, and that star-formation rates are derived from 
the far-infrared luminosity, it is important to determine the 
far-infrared spectral energy distributions for individual ob- 
jects if we are to obtain independent measures of the star- 
formation properties and potential of the quasar host galax- 
ies. We have started followup programs at several wave- 
lengths: (a) a parallel program at 1.25mm using the IRAM- 
30m which, when combined with (b) 450/^m measurements 
using SCUBA at the JCMT will be used to constrain the 
spectral energy distributions and thus dust temperatures 
and far-infrared luminosities (c) CO observations using mil- 
limetre interferometers to measure the gas mass, simulta- 
neously setting upper limits to the physical extents of the 
high-z host galaxies, are also in progress along with (d) VLA 
observations to sear ch for the signatures of starformation or 
weak radio-sources (Carilli et al. 200TI). 

Our approach to targeting starformation at high- 
redshift has been to preselect candidates that are likely to 
have massive molecular gas reservoirs, as inferred from the 
presence of dust. At some redshift, however, it may become 
more efficient to search directly for the gas reservoir itself, 
as the presence of significant quantities of dust preselects 
for objects in which at least one massive starburst of star- 
formation must already have taken place. 

A second selection effect arises from the original def- 
inition of our sample. By basing our sample on the most 
optically-bright quasars, we tacitly assume that the quasar 
and the host galaxy are co eval . Using arguments based on 

above, we have shown that 



timescales given in Section 6.3 



expect for an Arp220-like ultraluminous infrared galaxy at 
z ~4.5 (Figure ^ - the value obtained is between the curves 
predicted by the two different cosmological models). If we 
assume a source density for z > 4 quasars of 1 per 500deg -2 , 



this may be the case. The evidence, however, is not conclu- 
sive, and indeed depends on a number of assumptions, not 
least of which being that the rest-frame far-infrared lumi- 
nosity is powered solely by star- formation, and that there is 
no gravitational lensing of the quasar and its host galaxy. 
With the combination of measurements of the CO, dust and 
optical properties, the physical parameters that we can de- 
rive from these measurements and a sample of statistically 
significant size, we will be able to investigate in a more quan- 
titative manner. 

We are currently undertaking a parallel search for 
submillimetre-bright quasar host galaxies at redshifts z ~ 
2. The elapsed time interval between z ~ 4 and z ~ 2 is 
around 1.5 Gyr, corresponding to a doubling of the age of 
the Universe at z ~ 4. Given the gas masses observed in 
these radio-quiet quasar host galaxies, and the high star- 
formation rates that we infer from the FIR luminosity (Ta- 
ble 3), one might expect that if the formation redshift of 
quasar host galaxies is high, a sizeable fraction of any molec- 
ular gas reservoir would be used up between z <; 4 and by 
z ~ 2. The submillimetre-wave properties of the two sam- 
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pies should therefore be quite different. We note, however, 
that the e-folding lifetime of the quasar (50Myr) is consid- 
erably shorter than this interval (again, see Table 3), and 
so the quasar population that we see at z ~ 2 may not 
be the same as the z ~ 4 population. Any global evolu- 
tion in submillimetre-wave luminosity between z ii 4 and 
z ~ 2 host galaxies could in this case reflect variation in 
the underlying cosmological processes, such as merger rates, 
dynamical timescale and ava ilability of cold gas, driving t he 
evolution of the population (Kauffmannfc Haehnelt 2000). 
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8 SUMMARY 

We have observed a total of 38 zj;4 radio-quiet quasars at 
850/im using the JCMT. We have identified a total of 8 new 
submillimetre-bright sources with Sssopm > 10 m Jy includ- 
ing a previously ambiguous detection made at the IRAM- 
30m at 1.25 mm. The observations were carried out under 
modest observing conditions (sky transmission ~ 65%) , but 
even so the integration time per target was around 0.5 hour, 
giving an overall detection efficiency of 0.5 sources per hour 
of integration. This method, observing preselected targets in 
photometry mode, thus compares favourably with the com- 
plementary strategy of blank-field mapping. These observa- 
tions triple the number of known z<;4 objects for which CO 
observations should be possible in a finite amount of observ- 
ing time to N > 10. As such, we now have a submillimetre- 
bright sample of statistically significant size, for which we 
possess accurate identifications and spectroscopic redshifts, 
for which we can now start to study the star-formation prop- 
erties in detail. We analysed a new optical spectrum of the 
BAL PSS J1048+4407, one of the sources detected, and find 
that the redshift is quite differe nt fro m that used for CO 
observations of Guilloteau et al.,( 1999), which may explain 



the observed lack of CO from this object. The strong 850/ira 
emission in contrast to the lack of 1.35mm emission detected 
in PdB observations suggests strongly that the much higher 
angular resolution obtained with the interferometer is caus- 
ing significant fraction of the emission to be resolved out - 
the (sub)millimetre- source is extended. 

Of the remaining 30 objects observed, 27 were observed 
down to a sensitivity of 3<r < 12mJy and not detected. Based 
on this sample, we find that our detect ion rate is compara- 
ble to that measured by Omont et al.,( 1996b) at 1.25 mm 
using 30-m, scaled to match the appropriate flux limits. At 
present, the sample (detections and non-detections) is too 
small to determine formally whether there is a range of ab- 
solute magnitudes and redshifts for which the rate of detec- 
tion is high, however there may be some indication that the 
magnitude distribution of the detections and non-detections 
differ. We determined a fiducial flux of 2.0 ± 0.6mJy for the 
non-detections, which we believe to be a statistical mea- 
sure of the submillimetre flux of the more submillimetre- 
faint quasar host population. 
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